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Abstract 

We perform return interval analysis of 1 -min realized volatility defined by the sum of absolute high-frequency intraday 
returns for the Shanghai Stock Exchange Composite Index (SSEC) and 22 constituent stocks of SSEC. The scaling 
behavior and memory effect of the return intervals between successive realized volatilities above a certain threshold 
q are carefully investigated. In comparison with the volatility defined by the closest tick prices to the minute marks, 
the return interval distribution for the realized volatility shows a better scaling behavior since 20 stocks (out of 22 
stocks) and the SSEC pass the Kolmogorov-Smirnov (KS) test and exhibit scaling behaviors, among which the scaling 
function for 8 stocks could be approximated well by a stretched exponential distribution revealed by the KS goodness- 
of-fit test under the significance level of 5%. The improved scaling behavior is further confirmed by the relation 
between the fitted exponent y and the threshold q. In addition, the similarity of the return interval distributions for 
different stocks is also observed for the realized volatility. The investigation of the conditional probability distribution 
and the detrended fluctuation analysis (DFA) show that both short-term and long-term memory exists in the return 
intervals of realized volatility. 
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1. Introduction 

The study of extreme events has drawn much attention of scientists, for instance the nature records of floods, 
temperatures and earthquakes Jll 0. S B S Si ■ By investigating the return intervals between successive extreme 
events exceeding a certain threshold q, scaling behaviors are revealed in the return interval distributions for numerous 
complex systems QHSEl- As with this scaling behavior, we can infer that the probability distribution of the return 
intervals of rare events that are difficult to take empirical measurements. This scaling property of extreme events is 
supposed to be of great importance for the risk assessment of rare events. Further studies have shown that the scaling 
behavior of return intervals may arise from the long-term memory of the original records 13> H B [13 111- This 



suggests that the scaling behavior might also appear in other types of records with long-term correlations, such as the 
stock market records. 

The early study of extreme events in stock markets mainly concerns the return intervals between successive volatil- 
ities above a certain threshold q. The time intervals between consecutive trades and orders have also been widely 
studied lll^[T3i[T4ilT5i[l6ilT7ll . In general, the volatility R(t,6t) is simply defined as the magnitude of logarithmic 
return between tick prices at time t and t - 6t, i.e. R{t,6t) = | ^(^(f)) - ln(F(f - 6t))\, where Y{t) is the tick price 
at time f . With this definition, Yamasaki et al and Wang et al used the daily data and intraday data of US stocks to 
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study the probability distribution of volatility return intervals, and indeed found a scaling behavior |18, 2^ 21]. 
Further studies show that long-term memory also exists in the volatility return intervals. Similar scaUng behavior and 
long-term memory are observed in thousands Japanese stocks and 4 Chinese stocks ll22ll23[l . On the contrary, there is 
growing evidence showing that the return interval distribution may exhibit multiscaling behavior. Lee et al studied the 
1-min volatility data of the Korean KOSPI index ll24ll . and Wang et al analyzed the trade-by-trade data of 500 stocks 
composing the S&P 500 index and 1137 US common stocks |25, 2^, and Kaizoji analyzed the daily data of 800 
companies listed on Tokyo Stock Exchange and the Nikkei 225 index 12711 . and Ren and Zhou studied two Chinese 
stock market indexes |2^, and all these studies show that the return interval distributions for different thresholds q 
exhibit a systematic deviation from scaling and show multiscaling behavior Ren, Guo and Zhou conducted a more 
careful study to scrutinize 30 very liquid stocks in the Chinese market using the Kolmogorov-Smirnov (KS) test, and 
found that some stocks pass the KS test displaying scaling behaviors while others show multiscaling behaviors ll29ll . 

In the finance literature, there are many different estimators for volatility. Anderson et al proposed a daily realized 
volatility constructed from the sum of the square intraday returns t30Li31il . This realized daily volatility contains more 
valuable information about the intraday data, and is made arbitrarily close to the underlyin g in tegrated volatiUty. In 
addition, there are other daily volatility estimators constructed based on the intraday data 1321 133[ l34l l35n . Subse- 
quently, this realized volatility is generalized to the volatility which sums the square returns in a fixed time interval 
13611 . This realized volatility may better describe the price fluctuation caused by the trades occurring in that time 
interval. It is interesting to investigate the statistical properties of return intervals of realized volatility. 

In this paper, we investigate the statistical properties of the return intervals of realized volatility based on high- 
frequency intraday data in the Chinese stock market. Inspired by the sum of squared returns realized volatility raised 
by Anderson, we introduce an estimate of 1-min realized volatility by summing the absolute logarithmic returns 
utilizing all trading data in each minute. This 1-min realized volatility contains more valuable information which 
lies in the high frequency trading data. Using the Kolmogorov-Smirnov (KS) test and detrended fluctuation analysis 
(DFA) method, we test if the scaling behavior and long-term memory of the return intervals maintain with this 1-min 
realized volatility. 

The paper proceeds as follows: In Section |2] we introduce the database analyzed and the definition of 1-min 
realized volatility. In Section|3] we study the return interval distribution of the realized volatility using the KS tests. 
In SectionlH we further study the memory effect of the realized volatility return intervals. Section|5]summarizes. 



2. Volatility definition 



Our analysis is based on the high-frequency intraday data of SSEC and 22 liquid stocks traded on the Shanghai 
Stock Exchange. These 22 stocks are the most actively traded stocks representative of a variety of industry sectors, 
and consequently have the largest sizes among all the stocks. The prices and the associated times of the SSEC index 
and the individual stocks are recorded every six to eight seconds from January 2004 to June 2006. 

In many previous studies, the volatility is defined as the magnitude of the logarithmic return. 



7;i(f) = |ln}'(0-lnF(f-l)|, 



(1) 



where Y(t) is the closest tick price to a minute mark t EH I24I Eii I29I1 . In this paper, we focus on the realized 

volatility constructed from the sum of absolute trade-by-trade returns within one minute. Suppose that Y{t') is the tick 
price at time f', then the realized volatility is defined as 



^2(0= Yj llni'CO-lnFCf'-DI, 



(2) 



where the sum is taken over all the tick times between f - 1 and t. For both volatility definitions the sampling time is 
one minute, and the volatility data size is about 140,000 for each stock. 

Before doing the analysis, we removed the intraday pattern to eliminate its periodic effect on the return interval 
distribution 1 37 , 2^ via dividing the volatility /?,(f) by its average value corresponding to time t on the day | 2^ 21 ] . 
Then we normalize the data by dividing its standard deviation so that the volatility is in units of its standard deviation. 
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3. Probability distribution of realized volatility return intervals 



3.1. Empirical return interval distribution 

We study the return intervals t between successive volatilities exceeding a certain threshold q. For each threshold 
q, a series of return intervals is obtained and its empirical probability distribution Pqif) can be obtained. Many 
empirical studies have showed that the probability distribution function (PDF) of the scaled return intervals may obey 
a scaling form 

(3) 



where (t) is the mean return interval which depends on the threshold q. 

Ren et al have studied the return intervals of the volatility defined by R\ for the SSEC and the 22 constituent stocks 
analyzed in this paper |28, 3. They found that, for some of the stocks and the SSEC, the return interval distributions 
do not show scaling behavior In Fig. [i|(a) and (b), the PDFs of the scaled return intervals of volatility R\ for the 
SSEC and a representative stock 600028 are plotted. For the SSEC, the curves for different thresholds ^ = 2, 3, 4, 5 do 
not collapse to a single curve and show systematic deviation from scaling. For stock 600028, the deviation becomes 
relatively small, but one can still see some difference between the curves. Fig. [r|(c) and (d) plot the return interval 
PDFs of the realized volatility R2 for the SSEC and stock 600028 for comparison. We find that the curves for different 
q values approximately collapse to a single solid curve. This strongly suggests that the return interval distribution for 
the realized volatility shows better scaling behavior than that of the volatility defined by R\. 
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Figure 1: (Color online) Probability distributions of the scaled return intervals of volatilities defined by (a) R[ for SSEC, (b) R[ for stock 600028, 
(c) R2 for SSEC, and (d) i?2 for stock 600028. Solid curves in (c) and (d) are stretched exponential fits with y = 0.26 and 0.31 respectively. 



3.2. Kolmogorov-Smimov test of scaling in Pq(T)for different q 

Eyeballing of the return interval distributions offers a qualitative way of distinguishing the scaling and nonscaling 
behavior A quantitative method, the Kolmogorov-Smirnov (KS) test, is further adopted to examine the possible 
collapse of the interval distributions for different thresholds. We use the KS test to examine two return interval 
distributions for q - 2 and 5, which seems to behave most differently among all the distributions. If the distributions 
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for q - 2 and 5 pass the test, we can conclude that all the return interval distributions for different q values collapse 
onto a single curve and consequently obey a scaling law. Suppose that F2 is the cumulative distribution function 
(CDF) of return intervals for q - 2 and is the CDF of return intervals for ^ = 5. We calculate the KS statistic by 
comparing the two CDFs in the overlapping region: 

KS ^m^x{\F2-Fi\) . (4) 

When the KS statistic is smaller than a critical value denoted by CV (i.e., KS < CV), we can conclude that the 
distribution for ^ = 2 is coincident with the distribution for q - 5, and the return interval distributions show scaling 
behavior. The critical value is CV = c„/ ■\/mn/(m + «), where m and « are the numbers of interval samples for q, and 
qj Hllll, and the threshold is Ca - 1.36 at the significance level of a = 5% HEl. 



Ren, Guo and Zhou have used the KS test to examine the scaling behavior of the return interval distribution for the 
volatility defined by Ri, and found that among the 22 chosen stocks, 1 1 of the individual stocks pass the KS test and 
show scaling behavior, while the remaining 1 1 stocks as well as the SSEC do not show scaling behavior as depicted 
in Table [T] We also use the KS test to study the scaling behavior of the return interval distribution for the realized 
volatility defined by 7?2- Remarkably, only 2 stocks fail in the test, and the remaining 20 stocks pass the test and 
show scaling behavior More interestingly, the SSEC also exhibits scaling behavior, as manifested by the collapse of 
the distributions in Fig. [T] This confirms that the scaling behavior of the return interval distribution is significantly 
improved when the realized volatility R2 is adopted. The improved scahng behavior of the realized volatility return 
interval distribution is consistent with the scaUng behavior of the return interval distribution for 4 Chinese stocks 
reported by Qiu et al 1 23 3- 



3.3. Function fonn ofPq{T) 

Ren, Guo and Zhou have further performed the KS goodness-of-fit test 142114311 to study the particular form of the 
scaling function using the volatility defined hy R\, and confirmed that the scaling function of the stocks which show 
good scaling behavior could be approximated by a stretched exponential form 

/(t/(t)) = ce-'"^/<^»'. (5) 

We adopt the same method to test the hypothesis that the empirical distributions for different q values for the realized 
volatility defined by R2 are coincident with a single stretched exponential fit. Following Ren et al ll29ll . we test 
the return interval distributions for two special q values q - 2 and 5 which behave most differently among all the 
distributions. Only the 20 stocks which show scaling behaviors and the SSEC are tested. If both return interval 
distributions for q - 2 and 5 are identical to a same stretched exponential in the overlapping region of the scaled 
return intervals, we can conclude that the scaling function has a stretched exponential form. 

In the case of the KS goodness-of-fit test, the KS statistic calculates the difference between the cumulative distri- 
bution Fq of empirical return intervals and the cumulative distribution Fse from the fitted stretched exponential, 

KS = max [\Fq - FseI) ,^ = 2,5. (6) 



A weighted KS statistic, which is more sensitive on the edges of the cumulative distribution, is defined as 143 1 



KSW^mJ^^=I^\. (7) 

We generate 1000 synthetic samples from the best fitting distribution, and calculate the KS and KSW statistics for 
the synthetic data by taking the same measurements as we do for the empirical data as 

KSs.\m = max - /^sim.SEl) (8) 



'We have discussed with Qiu and revealed that they defined the voIatiUty in a similar way as we did for the realized volatility R2. However, they 
did not clearly clarify this in their paper (23ll . 
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Table 1 : Kolmogorov-Smimov test of the return interval distributions by compaiing the statistic KS with the critical value CV at the 5% significance 
level for SSEC and 22 constituent stocks. 



Stock code 




volatility K 


1 




Volatility R2 




Ao 




Scaling? 


AO 


CV Scaling? 




All nr\ 
U.i i /U 


U.U4j4 


rso 


rv rv/i cc 


0.0502 


Yes 


oUUUUU 




rv rvcTQ 

U.Ujjo 


JNO 


rv rvO/1 ^ 


0.0481 


Yes 


ouuuiy 




U.Uo // 


les 


rv rv/i r\n 
U.U4U / 


0.0701 


Yes 


600026 


0.0277 


0.0511 


Yes 


0.0362 


0.0516 


Yes 


600028 


0.0670 


0.0653 


No 


0.0374 


0.0652 


Yes 


600030 


0.0580 


0.0520 


No 


0.0409 


0.0492 


Yes 


600036 


0.0436 


0.0576 


Yes 


0.0466 


0.0535 


Yes 


600073 


0.0501 


0.0500 


No 


0.0345 


0.0452 


Yes 


600088 


0.0540 


0.0497 


No 


0.0250 


0.0456 


Yes 


600100 


0.0462 


0.0500 


Yes 


0.0239 


0.0461 


Yes 


600104 


0.0149 


0.0573 


Yes 


0.0365 


0.0510 


Yes 


600110 


0.0911 


0.0506 


No 


0.0958 


0.0457 


No 


600171 


0.0377 


0.0530 


Yes 


0.0204 


0.0482 


Yes 


600320 


0.0515 


0.0506 


No 


0.0305 


0.0469 


Yes 


600428 


0.0475 


0.0510 


Yes 


0.0345 


0.0501 


Yes 


600550 


0.0198 


0.0498 


Yes 


0.0198 


0.0467 


Yes 


600601 


0.0218 


0.0637 


Yes 


0.0460 


0.0569 


Yes 


600602 


0.0167 


0.0561 


Yes 


0.0316 


0.0513 


Yes 


600688 


0.0421 


0.0566 


Yes 


0.0391 


0.0531 


Yes 


600770 


0.0647 


0.0478 


No 


0.0284 


0.0442 


Yes 


600797 


0.1179 


0.0647 


No 


0.1220 


0.0568 


No 


600832 


0.0610 


0.0474 


No 


0.0203 


0.0469 


Yes 


600900 


0.0574 


0.0556 


No 


0.0300 


0.0550 


Yes 
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and 



max 



i.SEl 



(9) 



where Fsim is the CDF of each simulated synthetic sample and /^sim.SE is the CDF of its best fit obtained from integrat- 
ing the fitted stretched exponential. The p-value is defined as the frequency that KS sim > KS or KS Wsim > KS W, and 
it can be regarded as the probability that the empirical distribution is consistent with its best fit. The tests are carried 
out for SSEC and 20 constituent stocks, and the resultant /7-values are listed in Table|2] 



Table 2: KS and KS W goodness-of-fit tests of the scaling function form of return interval distributions for q = 2 and 5 by comparing empirical 
data with the best stretched exponential fit and synthetic data with the best stretched exponential fit. The stocks marked with * pass the test using 
the KS statistic, and the stocks marked with ** pass the test using both KS and KS W statistics under the significant level of 5%. 



Code 


q 


Pks 


Pksw 


Code 


q 


Pks 


Pksw 


Code 


q 


Pks 


Pksw 


SSEC* 


2 


0.321 


0.159 


600000** 


2 


0.338 


0.277 


600019* 


2 


0.074 


0.042 




5 


0.203 


0.002 




5 


0.434 


0.213 




5 


0.530 


0.324 


600028 


2 


0.011 


0.001 


600030 


2 


0.168 


0.094 


600026** 


2 


0.630 


0.406 




5 


0.712 


0.841 




5 


0.018 


0.005 




5 


0.288 


0.288 


600036** 


2 


0.707 


0.635 


600073** 


2 


0.417 


0.274 


600088** 


2 


0.443 


0.120 




5 


0.209 


0.205 




5 


0.411 


0.352 




5 


0.066 


0.056 


600104 


2 


0.013 


0.006 


600100** 


2 


0.379 


0.236 


600171 


2 


0.016 


0.001 




5 


0.045 


0.017 




5 


0.872 


0.665 




5 


0.005 





600320** 


2 


0.744 


0.526 


600428 


2 


0.010 


0.006 


600550* 


2 


0.070 


0.042 




5 


0.110 


0.051 




5 





0.001 




5 


0.564 


0.531 


600601** 


2 


0.493 


0.351 


600602* 


2 


0.522 


0.361 


600688 


2 


0.037 


0.001 




5 


0.417 


0.149 




5 


0.130 


0.039 




5 


0.007 


0.002 


600770* 


2 


0.680 


0.441 


600832* 


2 


0.532 


0.246 


600900* 


2 


0.307 


0.256 




5 


0.067 


0.015 




5 


0.122 


0.006 




5 


0.122 


0.005 



By checking the /j-values for both distributions for q - 2 and 5, we can test the null hypothesis that the empirical 
PDFs can be fitted well by a stretched exponential. Consider the significance level of 1%, if at least one /j-value for 
= 2 or 5 of an individual stock is less than 1%, then the null hypothesis is rejected. Table |2] depicts the /j-values 
for both distributions for q - 2 and 5 for the volatility defined by R2- According to Tabled 17 stocks as well as 
the SSEC pass the goodness-of-fit test using the KS statistic, and 12 stocks pass the goodness-of-fit test using the 
KSW statistic. For the volatility defined by only 8 stocks (out of the 20 constituent stocks examined by the KS 
goodness-of-fit test) pass the goodness-of-fit test using the KS statistic and 7 stocks pass the goodness-of-fit test using 
the KS W statistic P29!]- Consider the significance level of 5%, 14 stocks as well as the SSEC pass the goodness-of-fit 
test using the KS statistic, and 8 stocks pass the goodness-of-fit test using the KS W statistic for the volatility defined 
by R2- In contrast, for the volatility defined hy Ri, 6 stocks pass the goodness-of-fit test using the KS statistic and 
5 stocks pass the goodness-of-fit test using the KSW statistic 1I29I1 . We find that the KS and KSW statistics provide 
very similar results for the volatility defined by but the KS W statistic is more sensitive than the KS statistic for 
the volatility defined by R2. In principle, the /?- values of a stock are larger when the scaling of PDFs for different q is 
more significant. 



3.4. Dependence ofy on q 

It has been shown that for the stocks which show multiscaling behaviors, e.g. SSEC, the PDFs could also be well 
approximated by the stretched exponential distribution but with different exponent y for different threshold q 12811 . 
Therefore, we assume that the PDFs for the constituent stocks could be approximated by the stretched exponential 
function, whether they show good scaling behaviors or not. We fit the PDFs for different threshold q using a stretched 
exponential form 

PAt) = be-"^'l'^'^^\ (10) 
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where a, b and j are dependent of q if there is no scaling. We investigate the relationship between the exponent y and 
the threshold q to further study the tendency of return interval distribution with q. 

In Fig. |2]the exponents y are plotted as a function of the threshold q for SSEC and three representative stocks 
60000, 600028 and 600030. Though the curve for the volatility defined by R\ fluctuates a little which may due to the 
fitting errors caused by fluctuations, it shows an approximate decreasing tendency with an increase of q. Regardless of 
the fitting errors, the fitted exponent y mainly dominates the shape of the return interval distribution. The decreasing 
tendency of y further confirms our previous findings that the return interval distributions for different q values dis- 
tinctly differ from each other and show systematic deviation from scaling 1I29I1 . Whereas for the volatility defined by 
R2, y slightly fluctuates but stays relatively constant in comparison with that of the volatility defined hy R\. A similar 
phenomenon is observed for other constituent stocks. This provides additional evidence that the scaling behavior of 
the return interval distribution for the realized volatility defined by R2 is significantly improved. 




Figure 2: (Color online) Exponent y vs. threshold q for (a) SSEC, (b) stock 600000, (c) stock 600028, and (d) stock 600030. 

3. 5. Similarity of Pqir) for different stocks 

We also study the similarity of the return interval distributions for different stocks. We fix the threshold q - 2, 
and see how the return interval distributions behave for different stocks. In Fig. |3](a), the PDFs of the scaled return 
intervals of the volatility defined by Ri for SSEC and five representative stocks are plotted. It is seen that the curves 
for different stocks differ from each other and do not collapse to a single curve. We then use the volatility defined 
by Ro to investigate the PDFs of the scaled return intervals for the SSEC and the 5 representative stocks. As shown 
in Fig. [3](b), the curves for different stocks have very similar shapes and seem to collapse on a single curve. This 
indicates that the PDFs for different stocks may follow similar scaling function if we use the volatility defined by R2. 

4. Memory in the return intervals of realized volatility 

4.1. Short memory in realized volatility return intervals 

The memory effect is another important characteristic feature of the stock markets. Empirical study has revealed 
that the memory effect universally exists in the return intervals of various stock markets for the volatility defined by 
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Figure 3: (Color online) Probability distributions of the scaled return intervals for SSEC and five representative constituent stocks when q = 2 for 
volatilities defined by (a) Si and (b) R2- The soHd line in (b) is the stretched exponential fit with y = 0.26. 



Ri 2^ 22, 23, 221 ■ We try to test if the memory effect retains in the return intervals when we use the reahzed 



volatility defined by Ri- 

We first investigate the short-term memory of the return intervals by calculating the conditional PDF f ^(t|to), 
defined as the probability to find an interval t immediately after an interval tq. To get better statistics, we study 
f ^(t|to) for a bin of tq. The entire interval sequences are arranged in an ascending order and partitioned into four bins 
with equal size. ^'^(tIto) for tq in the smallest and biggest quarter bins for SSEC and stock 600028 with volatility 
defined by R\ are illustrated in Fig. |4](a) and (b). The curves intersperse with each other for small scaled return 
intervals, and it is hard to distinguish the curve for tq in the smallest subset from that for tq in the biggest subset. 
In comparison with that of the volatility defined by R\, /"^(tItq) for the same stock index and the individual stock 
for the volatility defined by R2 are plotted in Fig. |4](c) and (d). One observes that the scaling behavior of Pq{T\TQ) 
is significantly improved for the volatility defined by R2 especially for small scaled return intervals: the curves for 
all thresholds for tq in the smallest subset and biggest subset approximately collapse onto two separate solid curves. 
P^(t|to) shows bigger probabilities for tq in the smallest (biggest) subset when r/ (t) is small (big), and this indicates 
that small (big) intervals tq tend to be followed by small (big) intervals t. This can be regarded as a proof of the 
existence of short-term memory in realized volatility return intervals, while all the Pq{T\TQ) curves for the shuffled 
data collapse to a single exponential curve (not shown). 

4.2. Long memory in realized volatility and its return intervals 

To further investigate the long-term memory of the realized volatility and its return intervals we use the detrended 



fluctuation analysis (DFA) method 1144 I45L 1461 147L 14811 . known as a general method of examining the long-term 
correlation in time series analysis. The DFA method computes the average fluctuation F{1) of the cumulative series 
y{t) - Ylt'=\ ^('') of data series x(f') as 

J Ni I 

F{i)^-YjYp^^'^-y^^'^'^^^ (11) 

where Ni is the number of windows with fixed / data points, and y{t) is a local linear estimation for y{t) in a certain 
window /. It is expected that F{1) scales with I as 

F{T) ~ r, (12) 

The DFA method provides an accurate estimation of long-range correlation which do not depend on the length of the 
time series, and the scaling exponent a is supposed to be equal to the Hurst exponent when a <\ ll49ll . Generally, for 
a > 0.5 the time series are long-term correlated, and for a = 0.5 the time series are uncorrected. 

Fig. |5]presents the detrended fluctuation functions F(l) of the realized volatility for SSEC and stock 600028. A 
crossover behavior is observed: for small scales of /, F{1) obeys a power law with a relatively small exponent; while 
for large scales of /, F{1) obeys a power law with a relatively large exponent. The exponent a is estimated to be 
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Figure 4: (Color online) Conditional PDFs of the scaled return intervals in the smallest 1/4 subset (filled symbols) and the largest 1/4 subset (open 
symbols) for volatilities defined by (a) Ri for SSEC, (b) Ri for stock 600028, (c) R2 for SSEC, and (d) R2 for stock 600028. The Solid lines in (c) 
and (d) are the stretched exponential fits. 

0.68 + 0.01 and 0.86 ± 0.02 for SSEC and 0.69 + 0.01 and 0.98 ± 0.03 for stock 600028 respectively in small scale and 
large scale regions, apparently larger than 0.5. Similar results are observed for other constituent stocks. Therefore, we 
can conclude that the realized volatilities are long-term correlated. 

We also compute the detrended fluctuation functions F{1) of the return intervals of realized volatility as shown 
in Fig. |5] and observe similar crossover behavior. The estimation of the crossover point which separates the two 
power-law regions becomes an important task, since it essentially affects the determination of exponent a. A simple 
least squares estimation method is applied to determine the value of the threshold by minimizing the square distance 
between and its best power-law fits in small scale and large scale regions. In fact we have used the same least 
squares estimation method to find out the crossover point for of the realized volatility. In Fig. |5](a) and (b), the 
solid lines are power-law fits in small scale and large scale regions, respectively. Apparently the crossover tends to 
appear at smaller scales when the threshold q increases as quantitatively illustrated in Fig. |6](a). 

As with this least squares estimation method, we can further test the relation between the exponent a of the return 
intervals and the threshold q, and see how the long-term memory of the return intervals varies with the change of q. 
In Fig. |6](b) and (c), the exponent a for SSEC and stock 600028 are plotted as a function of the threshold q. The 
curves fluctuate a little for relatively large q due to the poor statistics of reduced interval samples. In general, the Hurst 
exponents in both regions for small scales and large scales show decreasing tendencies when q increases. Though the 
long-term correlation of the return intervals is weakened with the increase of q, the exponent a for all thresholds is 
apparently larger than 0.5. For the shuffled realized volatility data, the exponent a of the return intervals displays 
a value close to 0.5. This indicates that the long-term memory of the return intervals may arise from the long-term 
memory of original volatility records. 
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Figure 6: (Color online) (a) Crossover point as a function of q for SSEC (circles) and stock 600028 (diamonds). Exponent a of the realized 
volatility return intervals for large scales (diamonds) and small scales (circles) for (b) SSEC and (c) stock 600028. 



5. Conclusion 

In summary, we have studied the statistical properties of the return intervals of 1-min realized volatility based 
on the high-frequency intraday data for the SSEC index and 22 liquid constituent stocks. The Kolmogorov-Smirnov 
test shows that 20 stocks (out of 22 stocks) and the SSEC exhibit scaling behaviors. We found that the scaling 
behavior of the return interval distribution of the realized volatility is significantly improved compared with that of 
the ordinary volatility defined by the closest tick prices to the minute marks. We further adopted the KS goodness- 
of-fit test using the KS and weighted KS W statistics to study the particular form of the scaling distribution, and 
found the scaling function for 8 constituent stocks can be well-approximated by a stretched exponential distribution 
/(t/{t)) — ce""*'^^^'^^''. We calculated the relation between the exponent y estimated from the stretched exponential fit 
of Pq{T) and the threshold q, and further demonstrated the improved scaling behavior of the realized volatility. The 
similarity of Pqir) for different stocks is also observed for the realized volatility. 

We then investigated the memory effect of the realized volatility return intervals for the SSEC and 22 constituent 
stocks. Short-term memory is revealed by the observation of the conditional probability distribution P^(t|to) which 
also shows good scaling behavior for the realized volatility. Using the DFA method, we found that long-term memory 
exists in both realized volatility and its return intervals, and the exponent a of the realized volatility return intervals 
shows a decreasing tendency with the increase of the threshold q. 
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